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1. INTRODUCTION 


Finding the zeros of a nonlinear equation, f(x) = 0, is a classical problem of numer- 
ical analysis. Analytic methods for solving such equations rarely exit, and therefore, 
one can hope to obtain only approximate solutions by relying on iteration methods. 
For a survey of the most important algorithms, some excellent textbooks are available, 
see [4,8,10]. The classical Newtons method 
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n=0,1,2,.... (1.1) 


Being quadratically convergent, Newton’s method is probably the best known and most 
widely used algorithm. Time to time the method has been derived and modified in a 
variety of ways. One such method derived from Newton’s method by approximating the 
derivative with non-derivative term of difference quotient is Steffensen’s method [9,11]. 
The method requires two evaluations of function and is quadratically convergent. The 
interesting iterative scheme is Steffensen’s method that has the following form: 


f^ (z«) 
In order to control the approximation of the derivative and the stability of the iteration, 


a Steffensen's type method has been proposed in [2], this approach is based on a better 
approximation to the derivative f'(z4) in each iteration. It has the following form: 


(flan + o | f )| f (28)) — f(8))/ on f (En) f En) 


After that, the paper [1] has extended the above result on Banach spaces, obtained 
its local and semi-local convergence theorems, and made its applications on boundary- 
value problems by multiple shooting methods. 
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Tn+1 = Tn 
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A family of fourth order methods free from any derivative, satisfying the highest con- 
vergence order were established in [12-14]. 


2. q- CALCULUS 


In the following, q is a positive number, 0 < q < 1. For n € N = {0,1,...}, 
k € Z* = (,2,...] anda, a1,..., ay € C, the q-shifted factorial, the multiple q-shifted 
factorial and the q—binomial coefficients are defined by 





n—1l k 
(aq)o:=1, (a:4)n :— [D [Q1 — a4). (a1, 02...) = T (oi. (2.1) 
j=0 j=l 
d 1 +1 
o 7 — pî ut 1 a ,0—" 
H Siad ial 0-«4)01-4 7) 0-4 ) (2.2) 
01q nig (d: q)n 
respectively. The limit, lim (a;q)n, is denoted by (a;q)oo. Moreover (a;q)n has the 


representation, cf. [5], 
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(a; a)n = JD [7] aha. (23) 
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k=0 
The g—Gamma function, [5,6], is defined by 








q:d c 
n) e aoo sec en (2.4) 
’ oo 
where we take the principal values of q? and (1 — q)!~*. In particular 
(q: a) 
Paral) = ü D" neN. 


Let u € C be fixed. A set A C C is called a ij—geometric set if for x € A, ux € A. 
Let f be a function defined on a q-geometric set A C C. The q-difference operator is 
defined by the formula 
f(x) — f(qz) 
D = cA 0}. 2.5 
da) = DIST, se A-} (2.5) 
If 0 € A, we say that f has q-derivative at zero if the limit 
^) — f(0 
im LET) — f(0) 
noo xq” 
exists and does not depend on x. We then denote this limit by D,f(0). The q- 


integration of F. H. Jackson [7] is defined for a function f defined on a q-geometric set 
A to be 





, TEA (2.6) 





b b a 
[ (ou =| Odat- | f(t) dat, a,b € A, (2.7) 
where Et 
/ F(t) dgt = Yo 0 -af(od), m6 A, (2.8) 
n=0 


provided that the series converges. A function f which is defined on a q-geometric set 
A, 0 € A, is said to be q-regular at zero if 


lim f(xq")- f(0), for every x € A. 
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The rule of q—integration by parts is 
f st Dui dae = (Fo)(a) ~ Jim ooa" f Dagle)fla2)dyx. 29) 


If f, g are q—regular at zero, the limp,_,..(fg)(aq") on the right hand side of (2.9) will 
be replaced by (fg)(0). The two variable polynomial qv; (x, a), x, a € C, are defined to 
be 
z"(a/z;q),, xz 0, 
polz,a):= 1, Ynlx,a):= (2.10) 


n(n—1) 


(-Dyur- z7 a”, uem. 





In [3], Annaby and Mansour gave q- Taylor series in the following forms 














ex Df) * " 
füxy- 2. T. a ph a) + T.) f Pn-1ı(x, qt) D; f (t) dat. (2.11) 
n=l k(k—1 DE flag") 
fle) = ep (a, 2) 
kO Pate) (2.12) 
LI. (x, qt)D^ f(t) dot 
+ T(n) LAE Pn-1\2, q q q's 
3. A q-STEFFENSEN-SECANT METHOD 
In the following we set en = £n — à, €; = yn — à, Zn = Zn + qf (Ln), 
Un = Xn — f(22)/ fn, Zn], where fla, b] = ae 
_ D,f(a) , a(1-a4)D2/(a) , (1 —¢)?(1+4)Di f(a) (3.1) 
E ae T4(3) | T,(4) l 
Dif(a@)  a(1—a)(2 * a)D;f(a) 
m nd e 
and 
— Dif(a) 
= T.) ; (3.3) 


Now, we state and prove our q-Steffensen-secant Theorem with fourth order conver- 
gence. 


Theorem 3.1. Let f : D — R be a real-valued function with a root a € D, D C R, and 
let xq be closed enough to a. If D (a); k = 1,2,3 exist, and D,(a) £0, then 


Fn; Yn] — fln; Yn] + Flen, tn] 


f? [£5.95] fn); ne N, (3.4) 





Tn+1 = Yn 
is fourth-order convergent, and satisfies the following error equation 


enti = A 1 B(14- 44) | A 1C(14- 44) — A? B(3 + 2qA 4-24? 22) | + O(e2), nen. 
(3.5) 
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Proof. Using the Taylor expansion in (2.11), we have 






































f(za) = 
Dg f(a) Dif (a 
po) em WE pe nea) (uni ad (3.6) 
3 Fla In 
= (£n — a)(za — qa) (£n — qa) + T. I 3(a, qt) D; f (t)dgt. 
Rearranging the above equation again gives: 
f(z4) = Aen + Be? + Ce? + O(e4), (3.7) 
that is 
(zn) = flan + af (@n)) = 
1 tn+qf (tn) D, f (a) | 
EA pala, qt) DIS (Odat ET (En — a + af en) 
D? f(a) | 
3f(a 
EQ (rn — at af) tn qu + af len) Nn — d'a afl) 
(3.8) 
H _Daf(a) 
= Oleh) + TET (en +a) 
D; f(a) 
+E Gy (en t af Gen + afm) "a — a) 
3f (a 
EP (en ars) + af len) +all — a) en + afen) Fa = 0) 
= A(en + af (2n)) + Blen + qf (@n))” + C(en + af (an))* + Olek). 
Thus, 
f(n) = 
A[1 + qAlen + B[1 + 3¢A + q?A7Je2+ (3.9) 
Cll + 44A + 3g? A? + d? A?] + 29B?[1 + A] e? + O(e]). 
Moreover, 
(En + af (n)) 7 f (zn) 
Flen, Tn] m 
qf (tn) (3.10) 





=A + B[2 + qA]es + |C[3 + 34A + q? 4?] 4 Ld e? + O(ei). 
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Therefore, 
o Fln 0 
g(&n) CE Tig ds] = 
O(ed) +e, — A 1 B[1 + qA]ez+ (3.11) 


A B?[1 + qA][2 + qA] - qA1 B? — A71C[2 + 3qA + q2 27] 3 








Consequently, 
Daf (a) La ) 
O (£n — à — g(zn)) 4 T,(3) (£n — a — g(za))(za — qa — g(za)) 
Dj f (a) 





" T.) (En — a — g(za))(zu — qa — g(za)) (za — Ga — 9(z«)) 


1 5 —9(tn) 4 
i y3(a, qt) D, f (t)dot 











J-L——— 
Pal Dafa) (3.12) 
= Oeh) + Gy (en IC) 
D; f (a) 
TG — g(@n))(€n RE qf (za) T a(1 m q))+ 
3 a 
PAPE en = aber) len — alan) + alt — des — as) + aC — 1) 
= A(en — 9(2n)) + B(en — g(tn))? + C(en — g(z«))? + Oleh). 
This means 
f (Yn) = Olen) + B[1 + aA]e; — 
|a + qA][2 + qA] — qB? — C2 + 33A 4 gan e. eM 
and 
e* = O(et) + A 1 B[1 + q4]e2- 
[ae + qA][2 + qA] —qA4 1B? —A"*Cl2 + 3qA + ean e? oe 
On the other hand 
Fp) — La) fan 
ed (3.15) 
=A + Ben + lc +A 1B7{1 + A] e? + O(e3). 
Hence 
Flen; Yn] = O(en)+ 
(3.16) 


A? + 2ABen + |2AC + B?[8 + ZI e + [anc +2A TBR + Al] e 
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But 
fn) = Fn) _ 








Fln, ya] = SF 
ented = Sen) + glen) bun 
A -- B(1-4- qA)en 4 lca h qA)? + A!B?(1--4gA4- 244?) e? + O(ei). 
So that 
f lyn, Tn] — Flen, Yn] + Flen, £n] 
H Yn = = 
( ) Plyn, £n] 
AT! + jaca + qA) — A™’B(3 + 2qA + ean e (3.18) 
| — 2A 3 BC(2 + qA) + A^ B?(5 +3qA + an) e? + Olet). 
If we multiply H (xn) by f(yn) we get 
A(2n)f(Yn) = H (£n) flyn, alen = 
h + jaca +qA) — A ?2B(34-24A + 244) e24- (3.19) 
| —2A?BC(2 + qA) + A 2 B? (5 + 39A + wan) ee + ole) gx 
Taking in consideration that x41 is nothing but y, — H(an)f(yn) we get 
In41 =Yn — H (£n) f (Yn) 
=n — Í + jaca +qA) — A ?2B(3 + 2qA + ean) eic (3.20) 
| 247860 + 44) + AB + 3A APP) ed + O(C4)] 
Thus 
en+1 = jaca + qA) — A7? B(3 + 2qA + 297A”) + Olen) e?e 
(3.21) 


= A !B[I + qA] jaca + qA) — A ?B(3-- 2qA + 27 A?) ef + O(e?). 











This completes the proof. 





In order to compare our new method with Steffensen’s method, we give the following 
example. 


Example: In this example we take 


f(x) = cos(x) — z. 
The root of f(x) is a = 0.7390851332. Then the sequence {£n }n 
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Tn+1 = Un 


gf us) )|( cost) (Yn) — rn) Eq(an) + idan) 
Ean) fs) |È eostn) — costen) ) Balan) f| 





af En) | ( costus) — an) Ets) a) 





| a(n) + Fæn)| [È coston) — costs) ) Ealen) af) a 


is fourth-order convergent, where 


Eq(%n) = cos (a cos(£n) + (1 — aen ) — (1 +q) cos(£n) + qz. , 


a( costen) - 2n ) ' 


Eq(Xn) 





Yn = Xn 


Taking zo = 0, for q = 0.5, we find 





X1 X2 X3 T4 
Our’s 0.8617217519 0.7399567610 0.7390851885 0.7390851332 
Steffensen’s 2.175342650 0.76343368 0.7613122807 0.7595358304 


Taking zo = 1.1, for q = 0.9, we find 














X1 T2 T3 T4 
Our’s 0.7063822168 0.7388491909 0.7390851206 .7390851333 
Steffensen’s 0.8296038833 0.8040964255 0.7902498570 0.7814206993 


Taking zo = 1.35, for q = 0.001, we find 














X1 X2 X3 T4 
Our’s 0.8144712303 0.74139713209 0.7390873914 0.7390851090 
Steffensen’s 0.7429374052  0.7428816874  0.7428275625 0.7427749629 
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